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ABSTRACT 

We report the discovery of WASP-21b, a new transiting exoplanet discovered by the Wide Angle Search for Planets (WASP) 
Consortium and established and characterized with the FIES, SOPHIE, CORALIE and HARPS fiber-fed echelle spectrographs. A 
4.3-d period, 1.1% transit depth and 3.4-h duration are derived for WASP-2Ib using SuperWASP-North and high precision photo- 
metric observations at the Liverpool Telescope. Simultaneous fitting to the photometric and radial velocity data with a Markov Chain 
Monte Carlo procedure leads to a planet in the mass regime of Saturn. With a radius of 1.07 Rj„p and mass of 0.30 Mj„p, WASP-2Ib 
has a density close to 0.24 py„,, corresponding to the distribution peak at low density of transiting gaseous giant planets. With a host 
star metallicity [Fe/H] of -0.46, WASP-2Ib strengthens the correlation between planetary density and host star metallicity for the five 
known Saturn-like transiting planets. Furthermore there are clear indications that WASP-21b is the first transiting planet belonging to 
the thick disc. 

Key words, planetary systems ~ techniques: photometric - Techniques : radial velocities 



1. Introduction 

Observations of planets that transit their host star represent cur- 
rently the best opportunity to test models of exoplanet structure 
and evolution. These last ten years, the photometric surveys have 
led to an increasing list of transiting planets. More than seventy 
transiting planetary systems have been identified from Super- 
Earth to Jupiter-like planets. The WASP project operates two 
identical instruments, at La Palma in the Northern hemisphere, 
and at Sutherland in South Africa in the Southern hemisphere 
and led recently to the detection of about 30% of known tran- 



* Based on observations made with the SuperWASP-North camera 
hosted by the Isaac Newton Group on La Palma, the FIES spectro- 
graph on the Nordic Optical Telescope, the CORALIE spectrograph on 
the 1.2-m Euler Swiss telescope on La Silla Observatory, the SOPHIE 
spectrograph on the 1.93-m telescope on Haute Provence Observatory 
and the HARPS spectrograph on the 3.6-m ESO telescope at La Silla 
Observatory under programs 081.C-0388, 082.C-0040, 084.C-0185. 
Tables of photometric data are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
[http://cdsw eb.u-strasbg.fr/cgi-bin/qc at ?J/A-l-A/| 



siting planets. Each WASP telescope has a field of view of just 
under 500 square degrees. The WASP survey is sensitive to plan- 
etary transit signatures in the light curves of stars in the mag- 
nitude range V ~ 9-13. A detailed description of the telescope 
hardware, observing strategy and pipeline data analysis is given 
in Pollacco et al. 02006 ). 

Here we present the WASP photometry of 
SWASPJ230958.23-I-182346.0, high-precision photometric 
follow-up observations with the RISE instrument on the 
Liverpool Telescope and high-precision radial velocity (RV) 
observations with the FIES (2.6-m NOT), SOPHIE (1.93-m 
OHP), CORALIE (1.2-m Euler) and HARPS (3.6-m ESO) 
fiber-fed echelle spectrographs. These observations lead to the 
discovery of WASP-21b, a hot-Saturn transiting exoplanet. 

2. Observations 

2.1. SuperWASP observation 

SuperWASP-North is a multi-camera telescope system located 
in La Palma and consisting of 8 Canon 200-mm f/1.8 lenses 
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Fig. 1. Un-binned SuperWasp folded light curve for WASP-21 
with the ephemeris P=4.32 days and T0=2454743.0419. 



each coupled to e2v 2048x2048 pixel back illuminated CCDs 
(Pollacco et al. I2006I I. This combination of lens and camera 
yields a field-of-view of 7". 8 x 7". 8 with an angular size of 
13.7 arcsec per pixel. The target SWASPJ230958.23+182346.0 
was monitored from 2006 July 22 to November 27, totaling 
3814 photometric measurements (available in electronic form at 
CDS). The pipeline-processed data were detrended and searched 
for transits with the methods described in Collier Cameron et al. 
(I2OO6I 1. yielding a detection of a periodic transit-like signature 
with a period of 4.32 days. Figure [T] shows the un-binned 
Super WASP phase-folded Ught curve of WASP-21 with the 
ephemeris P=4.32 days and T0=2454743.0419. 



2.2. Higher precision piiotometric follow-up 

On-ofF photometric observations were made with the 2-m 
Faulkes Telescope North (FTN) at Hawaii on 2th September 
2008 with the Pan-STARRS-z broad-band filter. These observa- 
tions allowed the transit signal to be confirmed. Follow-up pho- 
tometry was performed with the fast read-out camera RISE lo- 
cated on the 2.0-m Liverpool Telescope at La Palma. Although 
the night was non-photometric and the sky was variable, a par- 
tial transit was observed on 7th October 2008 in agreement with 
the SuperWasp and FTN photometry. The phase-folded light 
curve of FTN and RISE (available in electronic form at CDS) 
are shown in Fig.|2] 

On October 2009, out-of-transit high-angular resolution im- 
ages were provided on the 4. 2-m WHT with the Adaptive Optic 
NAOMI and the IR imager INGRID. Ks-band images have a 
full width at half maximum (FWHM) of 0.18 arcsec with circu- 
lar contours, excluding evidence of a visible companion on that 
spatial scale. The limit in magnitude is estimated to be Ks~13.2 
at 0.36 arcsec and to be Ks~16.8 away from the tail of the point 
spread function (PSF). 

2.3. Doppler follow-up 

Three spectra were first obtained at the end of December 
2007 with the FIES spectrograph at the 2.6-m Nordic Optical 
Telescope (NOT - La Palma). The combined spectrum indi- 
cated a dwarf star and the radial velocities showed no varia- 
tions greater than 100 ms"^ excluding a low-mass star com- 
panion. Radial velocity (RV) follow-up of WASP-21 was con- 



Fig. 2. Phase-folded light curve from the RISE (open circle) 
and FTN (black circle) photometry for WASP-21 superposed 
with the MCMC solutions. The RISE photometry was binned 
2-minutes in time. 



ducted with the echelle spectrograph SOPHIE at the 1 .93-m tele- 
scope of Observatoire de Haute Provence (Perruchot et al. 120081 
Bouchy et al. 2009), CORALIE at the 1.2-m Swiss Euler tele- 
scope of La Silla (Baranne et al., 1996 Queloz et al. 2000 Pepe 
et al. l2002] l. a nd HA RPS at the 3.6-m telescope of ESO La Silla 
(Mayor et al. 12003). A total of eleven SOPHIE measurements 
were made in July and August 2008 with about the same signal- 
to-noise ratio S/N=30 in order to minimize the charge transfer 
inefficiency (CTI) effect (Bouchy et al. 2009) although it is now 
corrected by the data reduction software (DRS). Some measure- 
ments were made with moon light but without significant effect 
on the RV measurements considering that the RV of the target 
(-89.4 kms"') was always far from the moon RV. The RV un- 
certainties were computed assuming photon noise plus 15 m s"' 
quadratic ally added in order to take into account the systematics 
of the high efficiency mode like the guiding and centering errors 
(Boisse et al. 2010 ) and wavelength calibration uncertainty. Five 
CORALIE measurements were made simultaneously with those 
of SOPHIE and eight additional ones were made one year later in 
June and September 2009. One CORALIE point was made close 
to the transit epoch (BJD=2454686.83) and may be aff'ected by 
the Rossiter-McLaughlin effect. We did not use it to determine 
the Keplerian solution. However, considering the small v sin / (< 
2 kms ') of the star, the expected amplitude of the RM effect 
should be lower than 20 ms \ which makes it difficult to de- 
tect with CORALIE or SOPHIE. Nine HARPS measurements 
were made in September and October 2008 to definitively secure 
the candidate considering the small RV amplitude observed by 
SOPHIE and CORALIE, and also to exclude blended eclipsing 
binaries scenario and to provide a high S/N spectrum for stellar 
parameters analysis (see Sect. 13.11 ). Six HARPS measurements 
were added in October 2009 to check for eccentricity and for 
additional companion in the system. 

The phase-folded RVs in Fig. [3] shows the SOPHIE (black 
circles), HARPS (red squares), CORALIE (green triangles) and 
FIES (blue open circles) points. All are phase matched pretty 
well with the photometry. The best Keplerian fit was obtained 
with an RV offset from the HARPS data of -Hi 1, -20.5 and +15 
ms ' for the CORALIE, SOPHIE and FIES data respectively. 
The semi-amplitude of the RV curve is 7^=37.2+1.1 ms"'. We 
adjust the Keplerian orbit fixing the eccentricity to zero. There 
is no obvious evidence for eccentricity, although it is difficult to 
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fit for the eccentricity with a small amount of data points from 
four different instruments because we already have three other 
free parameters (the zero-point offsets of the RV data sets). 

No significant RV drift was detected with HARPS during one 
year excluding additional Jupiter-like companion with a period 
of less than one year. 



Table 1. Radial velocities of WASP-21 



BJD 


RV 


±lcr 


-2400000 


(kms-') 


(kms-') 


FIES - 


2.6m NOT 




54462.336 


-89.504 


0.027 


54465.371 


-89.409 


0.036 


54466.364 


-89.444 


0.029 


SOPHIE 


- 1.93m OHP 


54659.5872 


-89.409 


0.018 


54662.5947 


-89.418 


0.018 


54663.5749 


-89.392 


0.017 


54665.5826 


-89.449 


0.019 


54666.5773 


-89.458 


0.017 


54667.5790 


-89.387 


0.017 


54668.5925 


-89.381 


0.020 


54679.5940 


-89.474 


0.018 


54683.6242 


-89.495 


0.018 


54685.5591 


-89.393 


0.018 


54688.5479 


-89.445 


0.018 


CORALIE- 1.2mEuler 


54624.9095 


-89.427 


0.013 


54661.7883 


-89.487 


0.012 


54662.8643 


-89.475 


0.014 


54684.8541 


-89.442 


0.012 


54686.8333 


-89.430 


0.014 


55001.9020 


-89.424 


0.014 


55009.8829 


-89.398 


0.015 


55013.8951 


-89.431 


0.015 


55093.6551 


-89.506 


0.016 


55095.6656 


-89.444 


0.013 


55096.6795 


-89.407 


0.014 


55098.6554 


-89.500 


0.013 


55130.6417 


-89.433 


0.015 


HARPS 


- 3.6m ESO 


54706.7464 


-89.4001 


0.0100 


54709.7655 


-89.4792 


0.0033 


54710.7472 


-89.4287 


0.0038 


54749.6611 


-89.4309 


0.0036 


54754.6434 


-89.4154 


0.0031 


54755.6635 


-89.4317 


0.0025 


54760.6340 


-89.4639 


0.0022 


54761.6280 


-89.4782 


0.0027 


54763.5963 


-89.4040 


0.0031 


55111.6377 


-89.4797 


0.0029 


55112.6342 


-89.4361 


0.0032 


55113.6229 


-89.4118 


0.0027 


55115.6429 


-89.4892 


0.0050 


55116.6385 


-89.4655 


0.0023 


55117.6204 


-89.4184 


0.0045 



Figure 4 shows the bisector span versus RV. The SOPHIE 
and CORALIE bisector spans have too large uncertainties to ex- 
clude a possible correlation with radial velocities. The HARPS 
bisector spans definitively show no correlation, and scenarios of 
blended eclipsing binaries are definitively invalided. 
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Fig. 3. Phase-folded radial velocities of WASP-21 obtained 
with SOPHIE (black cu'cles), HARPS (red squares), CORALIE 
(green triangles) and FIES (blue open circles) 
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Fig. 4. Bisector span versus radial velocities of WASP-21 ob- 
tained with SOPHIE (black circles), HARPS (red squares) and 
CORALIE (green triangles) without significant coiTelation and 
invalidating scenarios of blended eclipsing binaries. 



3. Results and parameters 

3.1. Stellar parameters 

A detailed spectroscopic analysis of the stellar atmospheric 
properties was made with the FIES and HARPS spectra. We 
merged the HARPS spectra into one high-quality spectrum in 
order to perform a detailed spectroscopic analysis of the stellar 
atmospheric properties. This merged spectrum was continuum- 
normalized with a low-order polynomial to retain the shape of 
the broadest spectral features. The total signal-to-noise ratio of 
the combined spectrum was about 100 per pixel. 

As previously done for our analysis of WASP-1 (Stempels 
et al. |2U07] ) and WASP-3 (Pollacco et al. 200 8), we employed 
the methodology of Valenti & Fischer ( I2005I I. using the same 
tools, techniques and model atmosphere grid. We used the pack- 
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age Spectroscopy Made Easy (SME) (Valenti & Piskunov[T996|l, 
which combines spectral synthesis with multidimensional 
minimization to determine which atmospheric parameters best 
reproduce the observed spectrum of WASP-21 (effective tem- 
perature Tgff, surface gravity log g, metallicity [M/H], projected 
vsinj, microturbulence vjnic and the macroturbulence vjnac). 
The parameters we obtained from this analysis are listed in 
Table |2] 

We interpolated the Girardi evolutionary tracks using a 12 
Gyr old star with re//=5800 and [M///]=-0.4 and found Mv = 
4.97. A reddening E(B - V) of 0. 1 was estimated from the equiv- 
alent width of interstellar NaD lines with the Munari & Zwitter 
d 19971 1 relationship. With an observed Vmag of 11.58, the dis- 
tance is estimated to be 242^3° pc. This value agrees with the 
distance of 235 pc derived through the angular diameter {6 - 
0.042 + 0.002 mas) estimated from the infrared flux method 
and the stellar radius from Table |3] It also agrees with the dis- 
tance of 230 + 30 pc based on the measured V band magnitude 
(y = 11.58 + 0.08) and the absolute V magnitude of a typical 
G3V type star (Grav il988l l. 

The properties of WASP-21 suggest it is a member of the 
thick disc component of the Galaxy. With the NOMAD cata- 
logue proper motions (Zacharias et al. I2004I I and the measured 
systemic radial velocity, we calculated its U, V, and W space 
motion to be -34.6, -66.0, +62.1 kms"' relative to the Sun, 
adopting a distance of 242 pc based on the measured V band 
magnitude (V - 1 1 .58 + 0.08) and the absolute V magnitude of 
a typical G3V type star (Gray ^19881). It is clear that WASP-21 
is lagging behind the Sun (which moves with the thin disc) with 
a high negative V velocity similar to the canonical thick disc 
value of ~ -45 kms ' (Momson et al. [T990] l. WASP-21 is also 
metal-poor with {MjH} - -0.4 + 0.1, again similar to the mean 
thick disc value of 0.6 dex. In addition, the abundance anal- 
ysis shows enhanced a/Fe ratios, which is also typical of a thick 
disc star (see Edvardsson et al. 119931 1. The alpha elements, Mg 
Ca, Ti, and Si have abundances of — 0.25, approximately twice 
the value of iron, {FejH} =-0.46, and the other iron peak ele- 
ments ([y///]=-0.32, [Cr///] =-0.44, [Mn///]=-0.75, [Co///]=- 
0.39 and [Nij H]--QA6). Furthermore, when we compare the ef- 
fective temperature {T^ff - 5800+100 K), metallicity and mean 
stellar density (p/po = 0.84) derived from the transit parameters 
of WASP-21 to theoretical stellar evolution models from Girardi 
et al. (2000) as described in Hebb et al. (|2009jl, we find the star 
to be evolved off the zero-age main sequence and similar in age 
to that of the thick disc (age ~ 12 + 5 Gyr). In order to quantify 
the probability that WASP-21 is a member of the thick disc, we 
compared its chemical and kinematic properties to the properties 
of similar model stars in the Besangon Galactic model (Robin 
et al. 2003 ) where the population class of the stars are known. 
We selected a set of model main sequence stars in a small vol- 
ume nearby the position of WASP-21 for comparison. We gen- 
erated more than 800,000 model stars in 55 realizations of the 
simulation with spectral types of F5V-K7V, distances of 150- 
600 pc, and positions within +5° in each direction of WASP-21 
(/ = 92,b - -38). The model stars come from all Galactic com- 
ponents and are not selected based on their velocities or metal- 
licities. Of the model stars with metallicities between -0.5 and 
-0.3 dex and within +15 kms"' of the calculated U,V, and W 
values of WASP-21, all have ages > 7 Gyr (population class 7 
and 8) which suggests WASP-21 is old. Furthermore, 92% of 
the matching model stars belong to the thick disc. In summary, 
WASP-21 has a high probability of being a thick disc member 
based on its kinematic properties, low metallicity, abundance 
patterns, and likely old age. 



Table 2. Stellar parameters for WASP-21 derived from FIES and 
HARPS spectroscopy. 



Parameter 


WASP-21 


GSC 


01715-00679 


2MASS 


J23095825+ 1823459 


RA (J2000) 


23:09:58.25 


DEC (J2000) 


+ 18:23:45.9 


V 


11.58+0.08 


Distance 


242^50 pc 


Teff 


5800 ± 100 K 




4.2 ± 0.1 




-0.4 + 0.1 


V sin I 


1.5 + 0.6km/s 




-89.45 kms-' 


Age 


12 + 5 Gyr 


log A(Li) 


2.19 + 0.09 


[NalH] 


-0.47 + 0.11 


[Mg/H] 


-0.28 + 0.08 


[Ca/H] 


-0.25 + 0.13 


[Ti/H] 


-0.28 + 0.10 


[Mn/H] 


-0.75 + 0.14 


[Fe/H] 


-0.46 ±0.11 


ISi/m 


-0.25 ± 0.12 


ISc/m 


-0.33 ±0.11 


wm 


-0.32 ± 0.10 


[Cr/H] 


-0.44 ±0.17 


ICo/H] 


-0.39 ±0.14 


[Ni/H] 


-0.46 ± 0.14 



3.2. Planet parameters 

To determine the planetary and orbital parameters the radial ve- 
locity measurements were combined with the photometry from 
WASP, RISE and FTN in a simultaneous fit with the Markov 
Chain Monte Carlo (MCMC) technique. The details of this pro- 
cess are described in Pollacco et al. (120081 ). Recent features and 
improvements were included. The linear de-trending of the op- 
tical light curves was made with respect to phase at each step in 
the MCMC chain. An initial fit showed that the orbital eccen- 
tricity was poorly constrained by the available data and nearly 
consistent with zero. We therefore fixed the eccentricity param- 
eter at zero. We derived the host star mass consistently within 
the MCMC code by comparing the mean stellar density (mea- 
sured from the shape of the transit) and the observed metallicity 
and effective temperature of the host star to an empirical relation 
defined by a set of detached eclipsing binaries with component 
masses and radii measured to high precision (< 3%) (Torres et al. 
2010). With this novel technique, we derived a mass for WASP- 
21 of M = 1.01 ± 0.03Mo. The mass and radius of WASP-21b, 
given in Table[3]with the other best-fit parameters, are then found 
to be Mp= 0.300+0.01 1 Mj^p and Rp = 1.07+0.06 Rjup, respec- 
tively. 

Finally, we note that the host star mass derived from the the- 
oretical stellar evolution models of Girardi et al. (120001 ) is signif- 
icantly less than what we find using the empirical relation, even 
though we used the same observed stellar density, temperature 
and metallicity. Upon further investigation, we found that the 
metallicity has a strong effect on the radii and effective tempera- 
tures of the theoretical stars, which is not present in the observed 
eclipsing binary systems. Specifically, a ~ 5800 K, solar metal- 
licity (Z = 0.019) theoretical star with an age of 5 Gyr has a pre- 
dicted mass of I.OM0, but at Z = 0.008 ([M/H] = -0.37), the 
predicted mass is ~ 10% less (O.9M0). However, the foureclips- 
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Table 3. WASP-21 system parameters and Icr error limits de- 
rived from the MCMC analysis. 



Parameter Value 



Transit epoch To [HJD] 


2454743.0419!"-""'' 


Orbital period P [days] 


4.322482 


+0.000019 
-0.000024 


Planet/star area ratio (Rp/R,)^ 


0.01082 


+0.00037 


Transit duration tr [days] 


0.1398 


-0.0040 


Impact parameter h [R,] 


23+" '2 

"■^-'-0.15 


Orbital inclination / [degrees] 


88.75 


+0.84 
-0.70 


Stellar reflex velocity A' [ m s"'] 


37.2 d 


: 1.1 


Orbital semimajor axis a [AU] 


"■"-^^-0.00044 


Stellar mass M, [Mq] 


1.01 ± 


0.03 


Stellar radius R, [R^] 


1.06 ± 


0.04 


Stellar surface gravity log g, [cgs] 


4.39 ± 


0.03 


Stellar density p, [po] 


0.84 ± 


0.09 


Planet mass Mp [Mj] 


0.300 ± 


0.011 


Planet radius Rp [Rj] 


1.07 ± 


0.06 


Planet density p,, [pj] 


0.24 ± 


0.05 



ing binary components in Torres et al. (120101) with temperatures 
within +100 K of WASP-21, have [Fe/H] measurements ranging 
from -0.1 to -1-0.24, but masses similar to within 5%. There is a 
clear difference between the predicted effects of metallicity on 
the properties of stars and the observed effects in eclipsing bi- 
naries. Using a stellar mass of 0.9Mq leads to a planetary mass 
and radius for WASP-21b of Mp= 0.29+0.01 Mjup and Rp = 
1.03+0.06 Rjup, respectively, ~ 3% less than our adopted value. 

4. Discussion and conclusion 

Figure |5] shows the mass-radius diagram of known gas-giant 
transiting exoplanets (with a mass of more than 0.1 Mjup). 
WASP-21b is one of the only five transiting gas-giant planets 
with a mass below 0.4 Mjup including HAT-P-12b (Hartman et 
al. 12009 1 . HD 149026b (Sato et al. 2005 ), CoRoT-8b (Borde et 
al.llpTO') and WASP-29b (HelHer et al. 2010). 

Figure |6] shows the density histogram of transiting gas-giant 
planets with a mass between 0.2 and 1.5 Jupiter masses (total- 
ing 50 exoplanets). The probability density function presents 
a clear asymmetric distribution peaking at the lowest densities 
(0.25 pjiip) and decreasing from 0.25 to 1.25 pjup. WASP-21b 
and HAT-P-12b, both with a density of ~ 0.24 pjup, are bloated 
planets with a density close to the mode of that distribution. This 
distribution is biased because the transit detection method (espe- 
cially from ground) favours planets with large radii, hence with 
lower density for a given mass. If large radii are predominantly 
the result of irradiation, the selection effects will favour stars 
with small a/Rt, because that gives a stronger irradiation of the 
planet. Small a/R^ also increases the transit detection probabil- 
ity. 

Only four transiting exoplanets with a mass in the range 0.2- 
1 .5 M jup have a density greater than or equal to Jupiter including 
the two Saturn-like planets CoRoT-8b and HD149026b as well 
as W ASP-7b (HelUer et al. 2008i and OGLE-TR-113b (Bouchy 
et al. 120041 ). all of them having a host star with metallicity 
[Fe/H] > 0. The density of the five Saturn-like transiting planets 
appears to be clearly coiTelated with the metallicity of their host 
stars ([Fe/H]^-0A6, -0.29, +0.11, +0.3 and +0.36 for WASP- 
21b andHAT-P-12b, WASP-29b, CoRoT-8b and HD 149026b re- 
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Fig. 5. Mass-radius diagram of transiting gaseous giant exoplan- 
ets (Mp > 0.1 Mjup). 



spectively), confirming and reinforcing the relation established 
by Guillot et al. (2006). If low host star metallicity also favours 
larger planetary radii, it may explain the difference observed be- 
tween the host star metallicity distribution from radial velocity 
surveys (50% with [Fe/H] > 0.1) and transit surveys (40% with 
[Fe/H] > 0.1). 

It is clear from previous studies that thick disc stars do 
form planets, which means that planet formation was going on 
when the Milky Way was still relatively young (i.e. as early as 
~ 12 Gyr ago) and that these planetary systems have survived the 
process (i.e. major merger with a satellite galaxy) that formed 
the thick disc. WASP-21b is the first transiting planet with a 
measured radius and definitive mass to orbit a member of this 
population. 

Reid et al. (120071 ) identify five previously known exoplanet 
host stars, which are likely thick disc members based on their 
subsolar metallicities, relatively large space motions (compared 
to the Sun) and high [a/Fe] ratios. More recently. Neves et al. 
(j2009) identified five planet-host stars from the 29 stars belong- 
ing to the thick disc from the HARPS GTO planet search pro- 
gramme indicating that approximately 17% of the thick disc 
stars are planet hosts. They found some indication that metal- 
poor ([Fe/H] < -0.2) planet-host stars originate preferentially in 
the thick disc, as also suggested by Haywood (2008). This in- 
dication is however based on small number statistics and was 
contested by a new spectroscopic analysis made by Gonzalez 
(2009), showing that according to the mass abundance of the re- 
fractory elements (Mg, Si, Fe), which is important for planet for- 
mation, thick disc and [Fe/H] -poor thin disc planetary host stars 
have similar distributions. This author recommends to use a so- 
called refactory index [Ref] rather than [Fe/H] to compare the 
statistics of planets around thin disc and thick disc stars. Because 
thick disk stars are alpha-enhanced, they might have less Fe than 
thin disk stars, but they have relatively more Mg and Si, so they 
are not as metal -poor in terms of their ability to form planets 
as they seem to be at first glance. For WASP-21b, this refractory 
index [Ref] is -0.35. 
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mass in the range 0.2-1.5 Mjup. 
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